"0","```r
# get site lat/longs
meta <- read.delim(file=\"data/pollen/All.site.data-withagemodel-finalv2.txt\", sep=\"\t\")
site.longlats <- meta[match(pollen.6000.sites, meta$Handle), c('Longitude', 'Latitude')]
# simple relationships with latitude and longitude
par(mfrow=c(2,2), mar=c(4,4,1,1)+0.1)
plot(H~site.longlats$Latitude, pch=16, xlab=\"Latitude\")
plot(S.obs~site.longlats$Latitude, pch=16, xlab=\"Latitude\")
```"
"0","```r
plot(Simp~site.longlats$Latitude, pch=16, xlab=\"Latitude\")
plot(H_even~site.longlats$Latitude, pch=16, xlab=\"Latitude\")
```"
"0","```r
par(mfrow=c(2,2), mar=c(4,4,1,1)+0.1)
```"
